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Abstract 



[ The present paper solves the problem of the group classification of the general Burgers' equation ut = f{x, u)u^ + g{x, u)uxx, 
, where / and g are arbitrary smooth functions of the variable x and u, by using Lie method. The paper is one of the few 
■ applications of an algebraic approach to the problem of group classification: the method of preliminary group classification. 

A number of new interesting nonlinear invariant models which have nontriviai invariance algebras are obtained. The result of 

the work is a wide class of equations summarized in table form. 
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1 Introduction 

It is well known that the symmetry group method plays an important role in the analysis of differential equations. 
The history of group classification methods goes back to Sophus Lie. The first paper on this subject is [1], where Lie 
proves that a linear two-dimensional second-order PDE may admit at most a three-parameter invariance group (apart 
from the trivial infinite parameter symmetry group, which is due to linearity). He computed the maximal invariance 
group of the one-dimensional heat conductivity equation and utilized this symmetry to construct its explicit solutions. 
Saying it the modern way, he performed symmetry reduction of the heat equation. Nowadays symmetry reduction 
is one of the most powerful tools for solving nonlinear partial differential equations (PDEs). Recently, there have 
been several generalizations of the classical Lie group method for symmetry reductions. Ovsiannikov [2] developed 
the method of partially invariant solutions. His approach is based on the concept of an equivalence group, which is 
a Lie transformation group acting in the extended space of independent variables, functions and their derivatives, 
and preserving the class of partial differential equations under study [3,4,5]. 

The investigation of the exact solutions plays an important role in the study of nonlinear physical systems. A 
' wealth of methods have been developed to find these exact physically significant solutions of a PDE though it is 
rather difficult. Some of the most important methods are the inverse scattering method [6], Darboux and Backlund 
transformations [7], Hirota bilinear method [7,8], Lie symmetry analysis [9,10,11], etc. The paper [12], based on the 
Lie group method, is investigated a very famous and important equation, which is the general Burgers' equation as 
the form 

Ut = aul + bu^x, (1.1) 

where u = u{x,t) is the unknown real function, a, 6 G M and ab =/= 0. In the present paper, we consider the general 
Burgers' equation as the form 



Ut = fix, u)ul + g{x, u)u, 



(1.2) 
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where u = u{x,t) is the unknown real function, / and g are arbitrary smooth functions of the variable x and u. 
Eq. (1.2) represents the Burgers' equation combining both dissipative and nonlinear effects, therefore appears in 
a wide variety of physical applications. So it is important to lucubrate the exact explicit solutions and similarity 
reductions for this equation [13,14]. Here, we get the preliminary group classification of Eq. (1.2) by means of 
Lie point symmetry, and the constructed optimal systems of subalgebras. The knowledge of the optimal system of 
subalgcbras gives the possibility of constructing the optimal system of solutions [2,15,16] and permits the generation 
of new solutions starting from invariant or non-invariant solutions. 



2 Symmetry Methods 

Let a partial differential equation contains p dependent variables and q independent variables. The one-parameter 
Lie group of transformations 

Xi = Xi + €£,\x,u)+0{e^); = + eip°'{x,u) + 0{e^), (2.3) 

dx' du 
where = -g-\^^Q, i = 1, . . . and = -^-^|^^g, a = 1,. . . ,q, are given. The action of the Lie group can be 

recovered from that of its associated infinitesimal generators, we consider general vector field 

^=E^^(-'-)|-+i:^"(-'-)a^- (2-4) 

i=l a=l 

on the space of independent and dependent variables. Therefore, the characteristic of the vector field X given by 
(2.4) is the function 

Q-(x,u^'^)=^^{x,u)-J2^\x,u)^, a=l,...,q. (2.5) 

The second prolongation of the infinitesimal operator 

d d d 

X = C\x,t,u)—+e{x,t,u)—+^{x,t,u)— . (2.6) 

obtained via the following prolongation formulas: 

d t d d d 



OUx OUt OUxt OUxx 

the coefficients are obtained by 

= D,Q + (,^Ux, + i''uu, ip'^ = D,{D,Q)+(,^Ux,j + ^''uuj, (2.7) 

where Q = Lp — S}ux — £,'^ut is the characteristic of the vector field X given by (2.5). For instance 

ip^ = Dx^ - UxDxS} - utDx^\ (2.8) 

= Dt^ - UxDte - utDt^, (2.9) 

<p^^ = £>,<^^ - UxxDx^^ - UxtDx^, (2.10) 

^-t = Dtip^ - UxxDte - UxtDtf, (2.11) 

where the operators Dx and Dt denote the total derivatives with respect to x and t: 

d d d d d d d d 

-L'a: = TT + U-r— + (/,,;,,. Tj h (/,,.( Tj h . . . , A = ^ + Uj— + Utt^ h UtxT. 

OX OU OUx OUt ot OU OUt OUx 
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By the theorem 6.5. in [17], X^'^^ut — f{x,u)u^ — g{x,u)uxx\\(^i 2) ~ ^" ^i^^e 

X'-^^Ut - f{x,u)ul - g{x,u)uxx] =>f* - (/xC^ + fu<f)ul - {9xi^ + 9uf)uxx - 2/¥''^ 
therefore we obtain the following determining function: 

- (/x?' + fu^K - {g.e + 5«¥')«xx - 2/(p^«, - 1 (^ 2) = 0. (2.12) 
In the case of arbitrary f{x,u) and g{x,u) it follows that 

^1 = ^ = = ^* = = 0, (2.13) 

or 

e = ^=o, e = c. (2.14) 

Therefore, for arbitrary f{x,u) and g{x,u) Eq. (1.1) admits the one-dimensional Lie algebra Qi, with the basis 
X. = |. (2.15) 

Qi is called the principle Lie algebra for Eq. (1.1). So, it is remained to specify the coefficient / and g such that 
Eq. (1.1) admits an extension of the principal algebra Qi. Usually, the group classification is obtained by inspecting 
the determining equation. But in our case the complete solution of the determining equation (2.12) is a wasteful 
venture. Therefore, we don't solve the determining equation but, instead we obtain a partial group classification of 
Eq. (1.1) via the so-called method of preliminary group classification. This method was applied when an equivalence 
group is generated by a finite-dimensional Lie algebra Qg. The essential part of the method is the classification of all 
nonsimilar subalgebras of . Actually, the application of the method is simple and effective when the classification 
is based on finite-dimensional equivalence algebra Qg. 

3 Equivalence transformations 

An equivalence transformation is a nondegenerate change of the variables t, x, u taking any equation of the form (1.1) 
into an equation of the same form, generally speaking, with different f{x,u) and g{x,u). The set of all equivalence 
transformations forms an equivalence group S'. We shall find a continuous subgroup of it making use of the 
infinitesimal method. 

We consider an operator of the group S'c in the form 

d d d d d 

Y = ^^{x,t,u)—+^'^{x,t,u)—+(p{x,t,u)—+n{x,t,u,f,g)—+v{x,t,u,f,g) — , (3.16) 

from the invariance conditions of Eq. (1.1) written as the system: 

ut- f{x,u)ul~ g{x,u)uxx = Q, (3.17) 
/t = 5t = 0, (3.18) 

where u and /, g are considered as differential variables: u on the space (x, t) and /, g on the extended space {x, t, u). 

The invariance conditions of the system (3.17) are 

Y^'^'^ {ut - f{x, u)ul - g{x, u)uxx) = 0, 

Y('^){f,)=Y^^){g,) = 0, (3.19) 
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where F^^^ is the prolongation of the operator (3.16): 

y(2) = y + a ^ 9 ^ ^td d ^ t^^ ^t^. (3_20) 

OUt OUxt OUxx oft ogt 

The coefficients (p^ , ip* , ip^*^ , ip^^ , (p^* are given in (2.7) and the other coefRcients of (3.20) are obtained by applying 
the prolongation procedure to differential variables / and g with independent variables {x,u). In view of (3.18), we 
have 

= Dti^l) - fMe) - fuDti^), = Dt{u) - gMe) - 9uDt{v>), (3-21) 

d 

where Df = —. So, we have the following prolongation formulas: 
at 

1^* = Ht- fx£.l - fuipt, V* = vt- gxil - gu'ft- (3.22) 

The invariance conditions (3.19) give rise to 

f/ = u* = 0, (3.23) 
that is hold for every / and g. Substituting (4.49) into (4.48), we obtain 

f^t=i^t= 0, ='Pt = 0- (3-24) 

Moreover with substituting (3.20) into (3.19) we obtain 

ip* - 2f{x, u)ux'f' - g{x, u)ip'''^ - nul - uuxx -iy = 0. (3.25) 

Substituting (3.27), (3.28) and (3.29) into invariance condition (3.25), and introducing the relation Ut = f{x,u)u% + 
g{x,u)uxx to eliminate ut we arc left with a polynomial equation involving the various derivatives of u{x,t) whose 
coefficients are certain derivatives of ^^,£,^,(p,i2, and u. We can equate the individual coefficients to zero, leading to 
the complete set of determining equations: 

(3.26) 
(3.27) 

(3.28) 
(3.29) 

2?^) - g^uu (3.30) 
so, we find that 

^^(a;) = a(.x), C^(t) = cit + C2, (p{x,u) = Ciu + b{x), 

M = -2/(ci - a(x)), jy ^ -g{ci - a'{x)), (3.31) 

with constants Ci,C2 and two arbitrary functions a{x) and b{x) such that b"(x) = Ci — a' (a;). 

We summarize: The class of Eq. (1.2) has an infinite continuous group of equivalence transformations generated by 
the following infinitesimal operators: 

d d d d d 

Y = a{x)— + (cit + C2)^ + {ciu + Ka;))— - 2/(ci - aix))— - g{ci - a'{x)) — , (3.32) 

Therefore the symmetry algebra of the Burgers' equation (1.2) is spanned by the vector fields 

|9 d ^ , d d d 

8 f) ri a ^^-^^^ 

n = aix)-+2faix)-+ga'ix)-, ^ = K^)^. 





=e{x) 




e 


=e{t) = 











V 


= -get + 




M 
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Table 1 

Commutation relations satisfied by infinitesimal generators in (4.35) 
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Table 2 

Adjoint relations satisfied by infinitesimal generators in (4.35) 










[,] 


Vi 


Y2 


F3 


F4 




F5 




Vi 


Y2 


F3 


F4 




F5 


^2 


Vi 


Y2 


F3 


F4 - S F2 




F5 


Yz 


Vi 


Y2 


Fa 


F4 - S F3 




Fb 


Yi 


Vi 


e'Y2 


e^Fs 


F4 




Fb 






Y2 


F3 


F4 




Fb 



Moreover, in the group of equivalence transformations there are included also discrete transformations, i.e., reflections 
t — > -t, X — > -X, u — > -u, f — > -/, g — > -g. (3.34) 



4 Preliminary group classification 

One can observe in many applications of group analysis that most of extensions of the principal Lie algebra admitted 
by the equation under consideration arc taken from the equivalence algebra gs- We call these extensions ^-extensions 
of the principal Lie algebra. The classification of all nonequivalent equations (with respect to a given equivalence 
group Gg,) admitting (^-extensions of the principal Lie algebra is called a preliminary group classification. Here, 
Gg is not necessarily the largest equivalence group but, it can be any subgroup of the group of all equivalence 
transformations. 

So, we can take any finite-dimensional subalgebra (desirable as large as possible) of an infinite-dimensional 
algebra with basis (3.33) and use it for a preliminary group classification. We select the subalgebra Q5 spanned on 
the following operators: 

d d d 

ox ot ou 

r,^ d d ^ , d d d d „„s 

''^ = 'dt^''d-u-'^df-'d-g^ = + '^df + 'd-g- ^'-''^ 

The communication relations between these vector fields is given in Table 1. To each s-parameter subgroup there 
corresponds a family of group invariant solutions. So, in general, it is quite impossible to determine all possible 
group-invariant solutions of a PDE. In order to minimize this search, it is useful to construct the optimal system of 
solutions. It is well known that the problem of the construction of the optimal system of solutions is equivalent to 
that of the construction of the optimal system of subalgebras [2,15]. Here, we will deal with the construction of the 
optimal system of subalgebras of 35. 

Let G be a Lie group, with g its Lie algebra. Each element T Q G yields inner automorphism Ta — > TTaT~^ of 
the group G. Every automorphism of the group G induces an automorphism of g. The set of all these automorphism 
is a Lie group called the adjoint group G^. The Lie algebra of G^ is the adjoint algebra of g, defined as follows. Let 
us have two infinitesimal generators X,Y G L. The linear mapping KdXiY) : Y — > [X, Y] is an automorphism of 
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g, called the inner derivation of q. The set of all inner derivations a.dX{Y){X, Y E q) together with the Lie bracket 
[AdX, AdY] = Ad[X, Y] is a Lie algebra called the adjoint algebra of g. Clearly g"* is the Lie algebra of G^. Two 
subalgebras in g are conjugate (or similar) if there is a transformation of G"^ which takes one subalgebra into the 
other. The collection of pairwise non-conjugate s-dimensional subalgebras is the optimal system of subalgebras of 
order s. The construction of the one-dimensional optimal system of subalgebras can be carried out by using a global 
matrix of the adjoint transformations as suggested by Ovsiannikov [2]. The latter problem, tends to determine a list 
(that is called an optimal system) of conjugacy inequivalent subalgebras with the property that any other subalgebra 
is equivalent to a unique member of the list under some element of the adjoint representation i.e. i) Ad(g) f) for some 
g of a considered Lie group. Thus we will deal with the construction of the optimal system of subalgebras of 05. 
The adjoint action is given by the Lie series 

Ad(exp(s Yi))Yj =Yj-s [Y^, Yj] + ^ [Yi, [Y^, Yj]] (4.36) 

where s is a parameter and i,j = 1, • • • ,5. The adjoint representations of 55 is listed in Table 2, it consists the 
separate adjoint actions of each element of 35 on all other elements. 

Theorem 4.1. An optimal system of one- dimensional Lie subalgebras of general Burgers' equation (1-2) is provided 
by those generated by 

-Yi +Y^ = -dt + a, + 2fdf + gdg, 
Y4^ + Y^= tdt +dx + udu, 
-Yi + Fg = -tdt + dx- udu + 4fdf + 2gdg, 
Yi+Yi + Y5 = it + l)dt + dx+ udu, 
-Yi + Yi + Y5 = {t- l)dt + dx + udu, 

- I4 + n = (1 - t)dt - udu + 2fdf + gdg, (4.37) 
-Fi - F4 + n = -(1 + t)dt + dx- udu + ^fdf + 2gdg, 



1) 




Yi = dt. 


11) 


yll 


2) 


y2 = 


Y2 = d„ 


12) 


yl2 


3) 


y3 = 


Ys = du. 


13) 


yl3 


4) 


Y^ = 


Y4 = tdt + udu - "^fdf - gdg. 


14) 


yl4 


5) 


Y^ = 


Y5=dx + 2fdf + gdg, 


15) 


yl5 


6) 


Y^ = 


Yi+Y2=dt + dx, 


16) 


yl6 


7) 


Y^ = 


-Yi +Y2 = -dt + dx. 


17) 




8) 


y8 = 


Yi+Yi = it + l)di + udu - 2fdf - 


gdg, 




9) 




-Yi +Yi = it- l)di + udu - 2/5/ 






LO) 




--Y,+Y5 = dt + dx + 2fdf+gdg, 







Proof. Let 04 is the symmetry algebra of Eq. (1.2) with adjoint representation determined in Table 2 and 

Y = aiFi + 02^2 + 03^3 + aiYi + agYg, (4.38) 

is a nonzero vector field of g. We will simplify as many of the coefficients a^; i — 1, . . . , 5, as possible through proper 
adjoint applications on Y. We follow our aim in the below easy cases: 
Case 1: 

At first, assume that 05 7^ 0. Scaling Y if necessary, we can assume that as = 1 and so we get 

Y = aiYi + 02^2 + a3>3 + a^Y^ + Y^. (4.39) 
Using the table of adjoint (Table 2) , if we act on Y with Ad(exp(a2l2)), the coefficient of Y2 can be vanished: 

F' = aiFi + 03^3 + 04^4 + 15- (4.40) 
Then we apply Ad(exp(a3l3)) on Y' to cancel the coefficient of Y^: 

y" ^aiFi +04^4 + ^5. (4.41) 
Case la: 

If ai, 04 7^ then we can make the coefficient of Yi and I4 either -|-1 or —1. Thus any one- dimensional subalgebra 
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generated by Y with 03, 04 7^ is equivalent to one generated by ±Yi ±14 + 15 which introduce parts I4), 15), 16) 
and 17) oi the theorem. 
Case lb: 

For ai = 0, 04 7^ we can see that each one-dimensional subalgebra generated by Y is equivalent to one generated 
by ±14 + Ys which introduce parts 12) and 13) of the theorem. 

Case Ic: 

For ai ^ 0, 04 = we can see that each one-dimensional subalgebra generated by Y is equivalent to one generated 
by zbYi + I5 which introduce parts 10) and 11) of the theorem. 
Case 2: 

The remaining one-dimensional subalgebras are spanned by vector fields of the form Y with as = 0. 
Case 2a: 

If a4 ^ then by scaling Y, we can assume that 04 = 1. Now by the action of Ad(exp 0212)) on Y, we can cancel 
the coefficient of I2: 

Y ^a^Yi+azY^ + Yi. (4.42) 
Then by applying Ad(exp(a3l3)) on Y the coefficient of I3 can be vanished and we have 

y' =arYr + Yi. (4.43) 

The one-dimensional subalgebra generated by Y is equivalent to one generated by ztYi + 14 which introduce parts 
8) and 9) of the theorem. 
Case 2b: 

Let 04 = then Y is in the form 

Y = oiFi + 02^2 + asls- (4.44) 
Suppose that 02 ^ then if necessary we can let it equal to 1 and we obtain 

Y' = aiYi + Y2 + a^Ys. (4.45) 

By acting Ad(exp(a3y3)) on Y' , it changed to oiYi -I- Y2 : 
Case 2b- 1: 

Let ai be nonzero. In this case we can make the coefficient of Yi in Y either +1 or —1 and find 6), 7) sections of 
the theorem. 

Case 2b-2: 

If ai is zero then I2 is remained. Hence this case suggests part 2). 
Case 2c: 

Finally if in the latter case 02 be zero, then no further simplification is possible and then Y is one of cases of 
(4.37). There is not any more possible case for studying and the proof is complete. □ 

The coefficients f,g of Eq. (1.2) depend on the variables x,u. Therefore, we take their optimal system's projections 
on the space {x,u,f,g). The nonzero in x-axis or u-axis projections of (4.37) are 



(4.46) 



Proposition 4.2. Let Qm '■= {Yi, . ■ . ,Ym), be an m- dimensional algebra. Denote by Y^{i = l,...,r, < r < 
m, r S N) an optimal system of one- dimensional subalgebras of Qm and by {i = 1, • • • ,t,0 < t < r,t G N) the 
projections ofY^, i.e., Z^ = pr(y'). If equations 

f = ^x,u), g = <i/{x,u), (4.47) 



1) z^ 


= y2 


= y6 


= Y^ = d,, 


4) Z^ 


= y5 = 


= yio = yii = g^ 


+ 2fdf + gdg, 


2) Z2 


= y3 


= du, 




5) Z^ 




= yi4 = yi5 = 


+ udu. 


3) Z3 


= y4 


= y8 


= Y^ = -yi6 = udu - 2fdf - gdg, 


6) Z^ 


= yi3 


= yi7 = d^- udu 


+ 4fdf + 2gdg. 
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Table 3 

The result of the classification 



N 


Z 


Invariant 


Equation 


Additional operator X^^' 


1 




u 


Ut = ^ul + ^Uxx 


dx,dt + dx, -dt +dx 


2 


Z^ 


X 


Ut = ^ul + 'i'Ua:x 


du 


3 


z" 


X 


Ut = U^^ul + U'i'Uxx 


tdt + udu, {t + l)dt + udu, {t — l)dt + udu 


4 


z^ 


u 


Ut = e'''<i'ui+e'''i'Uxx 


dx, dt + dx, -dt + dx 


5 




u 


Ut = + U^Uxx 


tdt + dx + udu, {t + l)dt +dx + udu 


6 




_1_ 

U 


Ut = e'^^^ul + e'^-^Uxx 


—tdt + dx — udu, —(1 + t)dt + dx — udu 



are invariant with respect to the optimal system then the equation 

Ut = u)ul + u)u^^, (4.48) 
admits the operators X* = projection of on {t, x, u) . 
Proposition 4.3. Let Eq. (4-48) and the equation 

Ut = ^'{x,u)ul + -^\x,u)u^^, (4.49) 

he constructed according to Proposition 4-2. via optim,al systems Z"^ and Z*', respectively. If the subalgebras spanned 
on the optimal systems Z'' and , respectively, are similar in Qm, then Eqs. (4-48) and (4-49) are equivalent with 
respect to the equivalence group Gm, generated by Qm- 

Now we apply Proposition 4.2. and Proposition 4.3. to the optimal system (4.46) and obtain all nonequivalent 
Eq. (1.2) admitting (f-extcnsions of the principal Lie algebra g, by one dimension, i.e.. equations of the form (1.2) 
such that they admit, together with the one basic operators (4.50) of g, also a second operator X'-^K For every case, 
when this extension occurs, we indicate the corresponding coefficients f,g and the additional operator X^'^\ 

We perform the algorithm passing from operators Z^ {i = 1, ■ ■ ■ ,6) to f, g and X^^^ via the following example. 
Let consider the vector field 

Z'^ = d^-udu + ^fdf + 2gdg, (4.50) 

then the characteristic equation corresponding to Z^ is 

du df dg 

'^ = -u = if = i^ ^'-''^ 

and can be taken in the form 

/l=«e^ l2 = 4^, h = -^. (4.52) 
Prom the invariance equations we can write 

h = $(/i), h = *(/i), (4.53) 
it follows that 

/ = e"'$(A), g = e"'$(A), (4.54) 
where A = 7i . 
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From Proposition 4.2. applied to the operator we obtain the additional operator X^'^^ 

-tdt+d^-udu, -{l+t)dt + da,-udu. (4.55) 

After similar calculations applied to all operators (4.46) we obtain the following result (Table 3) for the preliminary 
group classification of Eq. (1.2) admitting an extension fla of the principal Lie algebra fli. 

5 Conclusion 

In this paper, following the classical Lie method, the preliminary group classification for the class of general Burgers' 
equation (1.2) and investigated the algebraic structure of the symmetry groups for this equation, is obtained. The 
classification is obtained by constructing an optimal system with the aid of Propositions 4.2. and 4.3.. The result of 
the work is summarized in Tabic 3. Of course it is also possible to obtain the corresponding reduced equations for 
all the cases in the classification reported in Table 3. We omitted these for brevity. 
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